Theoretical stopping power of copper for protons using the shellwise local plasma 

approximation 
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We present a theoretical study on the energy loss by protons in solid copper. The ab initio 
shellwise local plasma approximation (SLPA) is employed for the inner-shells, and the Mermin 
dielectric function for the valence electrons. The partial contribution of each sub-shell of target 
electrons is calculated separately, including the screening among the electrons of the same binding 
energy. Present results are obtained using the SLPA and the known Hartree-Fock wave functions 
for copper. The theoretical curves are compared with the experimental data available and with 
the semi-empirical SRImOS values, showing very good agreement in a large energy region (5 keV-30 
MeV). 

PACS numbers: 



I. INTRODUCTION 

The stopping power of ions in solids is a necessary in- 
gredient of many areas of basic science and material tech- 
nology [U, 01 . There are numerous detailed models to de- 
scribe the response of outer electrons as a free electron 
gas (PEG) [1, 0, Q . In many cases the stopping due to 
the FEG is enough for energies reaching up to the max- 
imum of the stopping. But for targets with the full d or 
even / subshells, the bound electron contribution plays 
an important role even at intermediate energies. The 
description of bound electrons from first principles is a 
heavy task, especially for multielectronic atoms (see, for 
example the binary coUisional formalism in Refs. 

In this work we describe the bound electron contribu- 
tion to the stopping power by employing a many electron 
model, the shellwise local plasma approximation (SLPA). 
This is based on the local plasma approximation (LPA) 
d, H, to deal with bound electrons as an inhomoge- 
neous free electron gas. We apply this model within the 
dielectric formalism fll|, IT^ uM but assuming indepen- 
dent shell approximation [lj|. This means that in the 
SLPA the electrons of the same binding energy respond 
to the ion passage as a whole (collectively), screening the 
interaction with the impinging ion. This binding energy 
(or ionization gap) is included explicitly by employing 
Levine and Louie dielectric function [l5| . 

The SLPA allows us to calculate the different moments 
of the energy loss by the ion when at least one of the 
bound electrons is ionized. This is an ab initio calcu- 
lation - no parameter is included - whose only inputs 
are the atomic densities of the different subshells and the 
corresponding binding energies. Once the electronic wave 
functions are available, the SLPA applies to the different 
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targets with the same degree of complexity. This model 
has been applied successfully to describe experimental 
stopping cross sections of elements with Z < 54 (using 
the Hartree-Fock wave functions and binding energies) 
[iBl, and also for very heavy elements like Au, Pb, and 
Bi for which numerical solutions of the Dirac equation 
were needed [l7j . 

In what follows we present our theoretical results for 
the stopping of protons in Cu (Z=29, [Ar] 3d^°As^). The 
aim of this contribution is to use Cu as a testing case for 
the SLPA because it is a well known element (experimen- 
tally [li,!!^! and theoretically [13, HH, Ull) , and because 
we can employ the known Hartree-Fock wave functions 
[2^ . [23 | to obtain the subshell densities to include in the 
SLPA. 

We organize this work in two steps. First, we describe 
briefly the theoretical calculation, and then we present 
results for stopping cross sections of Cu for protons and 
compare our ab initio theoretical values with the exper- 
imental data available and with the semi-empirical 
SRIMO8 results il9|]- 



II. THEORETICAL CALCULATIONS 

The SLPA introduces two important differences to the 
known LPA. One of them is the shell-to-shell calcula- 
tion. Mathematically it implies a separate dielectric func- 
tion for each subshell of bound electrons, and the total 
response as the addition of these independent contribu- 
tions. Physically, the electrons are screened only by those 
of the same binding energy and not by the rest. 

The second one is the explicit inclusion of the bind- 
ing energy by using the Levine-Louie dielectric response 
[m instead of Lindhard one \2§\ . This dielectric function 
takes into account the ionization gap so that the contri- 
bution is null if the ionization threshold is deeper than 
the energy transferred The Levine-Louie dielectric 
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response maintains the characteristics of Lindhard [25| 
(Hnear response, electron-electron correlation), and satis- 
fies the /-sum rule (particle number conservation) . More 
detail about the SLPA can be foimd in Refs. [l6j. 

For a bare ion of charge Zp and velocity w, the stopping 
cross section due to the ionization of the nl subshell of 
target electrons, is expressed as 
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The dielectric function Sniik, to) is calculated as a mean 
value of a local response 
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where Sni{r) and i?„; 
(5„i(r) and ionization 



(2) 

are the local density of electrons 
gap, respectively, obtained from 
the atomic Hartree-Fock tabulated by Bunge et al. [1^ . 
The spatial integration is performed over the atomic di- 
mensions within the solid Cu (i.e., Rws is the atomic 
Wigner-Seitz radius, related to the atomic density 5at — 

The dielectric response, given by Eq. ([2]), considers 
only the electrons in the nl subshell. In this way the 
contribution of each shell of target electrons is obtained 
separately. An orbital version of the LPA was introduced 
many years ago by Meltzer et al. [26] within the logarith- 
mic high energy limit of the stopping number. Instead, in 
the SLPA we calculate a dielectric response of each shell, 
and it is valid even in the intermediate energy range. 

The independent shell description in the SLPA con- 
siders explicitly the energy threshold Eni of each shell 
(ionization of inner shells) . The Levine and Louie model 
[151] defines the dielectric function as 



The characteristic plasmon frequency and width of the 
Cu FEG are Wp = 0.703 a.u. and 7 — 0.950 a.u. respec- 
tively. This implies a mean value of electrons in the FEG 
Ne = 3.14, and a Seitz radius rs — 1.82 a.u. .These values 
were obtained from the optical data of energy loss func- 
tion (30j by considering only the first important peak. 
The number of electrons in the FEG is similar to the 
experimental value recommended by Isaacson [3l| . To 
keep the total number of electrons, we considered Cu as 
[Ar] 3d^ ®^ and 3.14 electrons as FEG. 
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Lindhard dielectric function [25|. This modified SLPA 
has already been applied successfully in recent calcula- 
tions of stopping power [ll, energy loss straggling 
in solids, and also in multiple ionization cross sec- 
tions of rare gases [14l] . 



FIG. 1; Stopping cross sections of Cu for protons. Curves: 
Solid-lines, present theoretical results for the contributions 
by bound electrons (SLPA) and the FEG, and total stopping 
as the addition of the previous two; dash-line, SRImOS results 
[l9l | . Symbols: black hollow- circles, sets of experimental data 
compiled by H. Paul [l^ corresponding to (1935-1979); dif- 
ferent symbols signed within the fi gure , experimental data in 

Refs. M H, Mm, M M M MM, M m, S, H, S, El . 



III. RESULTS 

In order to compare with the experimental data avail- 
able we calculate total stopping cross sections. These 
values are obtained as the addition of the bound elec- 
tron and the FEG contributions. The former calculated 
with the SLPA formalism as described before. The latter 
by employing the dielectric formalism with the Mermin 
dielectric function 12911. 



In Fig. [T] we plot our total stopping cross sections of 
Cu for protons, and compare them with the experimen- 
tal data available and with SRImOS results [ll|. The 
contribution from the FEG and the bound electrons are 
displayed separately. About the experimental data, we 
indicate separately only the data since 1980, while previ- 
ous one is plotted together with a single type of symbol. 
The agreement with the experimental data and with the 
SRImOS curve is good. This is a perturbative formalism. 
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which is not expected to describe the data for ion en- 
ergies below 30 keV. However, we have obtained a very 
good a gree ment even with the measurements by Valdes 
et al. 3z| in the low energy region. The FEG results 
obtained with the Mermin dielectric function, even per- 
turbative, are rather good for protons in Cu at such low 
energies. Figure [T] shows that the stopping maximum 
is correctly described in energy and value. For energies 
above that of the stopping maximum, the description is 
good, but shows an overestimation around 2 MeV. This 
behavior of the SLPA has already been noted in previous 
works [itI [13] and it is under study. 

IV. CONCLUDING REMARKS 

In this work we test the ab initio SLPA for the stop- 
ping of Cu for protons. This is a perturbative model 



which works within the dielectric formalism by consid- 
ering the response of bound electrons as that of a free 
electron gas of inhomogeneous density. The SLPA intro- 
duces the independent shell approximation to the usual 
formalism. The theoretical results are compared with the 
experimental data available and with the semi empirical 
SRImOS values, showing good agreement in a large energy 
region (5 keV-30 MeV). 
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